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ABSTRACT 

We calculate synchrotron radiation in three-dimensional pseudo-Newtonian magnetohydrodynamic 
simulations of radiatively inefficient accretion flows. We show that the emission is highly variable at 
optically thin frequencies, with order of magnitude variability on time-scales as short as the orbital 
period near the last stable orbit; this emission is linearly polarized at the ~ 20 — 50% level due to the 
coherent toroidal magnetic field in the flow. At optically thick frequencies, both the variability amplitude 
and polarization fraction decrease significantly with decreasing photon frequency. We argue that these 
results are broadly consistent with the observed properties of Sgr A* at the Galactic Center, including 
the rapid infrared flaring. 

Subject Headings: Accretion, accretion disks - Galaxy: center 



1. INTRODUCTION 

The identification of magnetohydrodynamic (MHD) tur- 
bulence as a robust angular momentum transport mecha- 
nism has led to a new era in the study of accretion onto 
compact objects (see Balbus 2003 for a review). Numer- 
ical simulations can now demonstrate accretion without 
imposing an arbitrary anomalous viscosity. Although an- 
alytic steady state models provide an important frame- 
work for understanding the structure of accreting systems 
(e.g., Shakura and Sunyaev 1973; Narayan and Yi 1994; 
Blandford and Begelman 1999), a computational approach 
is required to capture the time-dependent, turbulent evo- 
lution of the flow. Simulations open up the time domain 
to theoretical study and should ultimately provide a more 
robust description of the observed features of accreting 
systems. Krolik and Hawley (2002) exploited this possi- 
bility by computing power spectra of quantities such as the 
accretion rate and Maxwell stress, both of which are pre- 
sumably related to the bolometric output of the accretion 
flow. Hawley and Balbus (2002) carried out similar calcu- 
lations, including estimating the characteristic frequency 
for synchrotron radiation in their simulations. Here we ex- 
tend these calculations by explicitly calculating the radia- 
tive output from MHD simulations; for reasons explained 
below, we focus on synchrotron radiation. Our goal is 
to make general connections between numerical simula- 
tions and observations of accreting systems - we do not 
believe that the simulations yet contain sufficient physics 
to warrant detailed quantitative comparisons. We focus 
on observables such as polarization and fractional variabil- 
ity that are not as sensitive to the absolute value of the 
predicted flux (which is difficult to calculate for several 
reasons explained in §2). 

We apply our results to observations of Sgr A*, the 
« 3.6 X IO'^Mq massive black hole (MBH) at the cen- 
ter of our galaxy (Schodel et al. 2002; Chez et al. 2003). 
Sgr A* is known to be remarkably faint: the canonical 



Bondi accretion rate estimate inferred using the observed 
gas density around the black hole predicts a bolometric 
luminosity L « O.lMc^ k, 10''^ ergs s~^, if radiation were 
produced with 10% efficiency (e.g., Baganoff et al. 2003). 
Instead, the observed luminosity is ~ 10^^ ergs s~^ (Melia 
and Falcke 2001). This disagreement favors radiatively in- 
efficient accretion flow (RIAF) models of Sgr A*, in which 
little of the gravitational binding energy of the inflowing 
matter is radiated away (e.g., Narayan et al. 1995; see 
Quataert 2003 for a review). Numerical simulations of 
RIAFs have shown that part of the reason Sgr A* is so 
faint is that the accretion rate is much less than the Bondi 
estimate (e.g.. Stone & Pringle 2001; Hawley & Balbus 
2002; Igumenshchev et al. 2003). This conclusion has 
been confirmed by observations of linear polarization in 
the mm emission from Sgr A* (Aitken et al. 2000; Bower 
et al. 2003; see, e.g., Agol 2000, Quataert & Gruzinov 
2000, and Meha, Liu, & Coker 2000 for how the observed 
polarization constrains accretion models). 

Sgr A* has been detected from the radio to the X- 
rays, with most of the bolometric power emerging at wave- 
lengths near 1 millimeter (see, e.g., Melia & Falcke 2001 
for a review of spectral models). The millimeter emission 
is significantly linearly polarized (~ 10%), consistent with 
synchrotron radiation (Aitken et al. 2000; Bower et al. 
2003), while lower frequency radio emission is circularly 
polarized at the k. 0.1 - 1% level (Bower et al. 1999). For 
the purposes of this paper, the most interesting aspect of 
the observations is the strong variability detected across 
the entire spectrum. X-ray flares with amplitudes from 
factors of a few up to ^ 50 have been detected on hour 
time-scales (Baganoff et al. 2001; Porquet et al. 2003). 
Flaring has also been seen in the infrared on time-scales 
as short as 10s of minutes (Genzel et al. 2003; Ghez et 
al. 2004). Order unity variability in the millimeter has 
been detected on longer time-scales (week to month; Zhao 
et al. 2003), with even lower amplitude, longer time-scale 
variations in the radio (Herrnstein et al. 2004). Occasion- 
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ally, the radio emission can vary more rapidly as well, by 
~ 25% over a few hours (Bower et al. 2002). In addition 
to variability in the total flux, the polarization in the ra- 
dio exhibits variability as well (Bower et al. 2002, 2005). 
The origin of these variations is not fully understood. An- 
alytic models generally interpret the X-ray and IR flares 
as emission due to transient particle acceleration close to 
the black hole (e.g., MarkofF et al. 2001; Yuan, Quataert, 
& Narayan 2003, 2004; hereafter YQN03 and YQN04), al- 
though alternative interpretations exist (e.g., Liu & Melia 
2002; Nayakshin et al. 2004). 

In this paper we compute synchrotron radiation in an 
MHD simulation of a RIAF and discuss the results in the 
context of observations of Sgr A*. We focus our analysis on 
synchrotron radiation because it probably dominates the 
bolometric output of RIAFs at very low accretion rates 
(see, e.g., Fig. 5 of YQN04), and because of its clear ap- 
plicability to Sgr A* (at least in the radio-IR). The syn- 
chrotron radiation is calculated "after the fact," using the 
results of existing numerical simulations that do not in- 
corporate radiation. This is self-consistent because energy 
loss via radiation is negligible and does not modify the 
basic dynamics or structure of the accretion flow. 

The plan for the rest of this paper is as follows. In 
the next section (§2) we describe the simulations and our 
method for calculating synchrotron radiation. We then de- 
scribe our results in §3. We begin by considering optically 
thin synchrotron emission (both total intensity and polar- 
ized), and then consider optical depth effects. Finally, in 
§4 we conclude and discuss our results in the context of 
observations of Sgr A*. Before describing our method in 
detail, it is worth noting some of the simplifications in our 
analysis. The simulations used in this paper arc Newto- 
nian and so do not account for relativistic effects near the 
black hole. In keeping with this simplified dynamics, we 
do not account for Dopplcr shifts or gravitational redshift 
when calculating the radiation from the flow. Lastly, the 
MHD simulations used here do not have detailed infor- 
mation about the electron distribution function, which we 
parameterize as a Maxwellian with a power-law tail. 

2. SYNCHROTRON RADIATION FROM MHD SIMULATIONS 

The speciflc simulation used in this paper is 'Model A' 
described in Igumenshchev et al. (2003). This is a Newto- 
nian MHD simulation using the Paczynski-Wiita potential. 
The domain of the simulation is a set of 5 nested Cartesian 
grids, each of which has 32 elements in the x and y direc- 
tions and 64 elements in the z direction. The smallest grid 
has elements that are 0.5 Schwarzschild radii (Rg) across, 
where Rg = 2GM/c?, and each successively larger grid 
has elements twice as large as the previous grid. All of the 
grids cover only two octants, both positive and negative 
z, but only positive x and y. The boundary conditions are 
periodic such that the x = plane is mapped onto the y = 
plane. At small radii the black hole is represented by an 
absorbing sphere and the central 2 R^ are not evolved. At 
the outer boundary, R = 256 J?s, outflow boundary condi- 
tions are applied. Matter is injected continuously near the 
outer boundary with toroidal magnetic fields, allowing a 
steady state to be reached in which inflow into the grid is 
balanced by outflow out of the simulation domain and into 
the "black hole." The flow has zero net poloidal magnetic 



fleld, and forms a geometrically thick convective accretion 
disk. There are no jets or significant mass outflow in the 
polar direction. 

Since the simulation variables are dimensionless, length- 
scales and time-scales can be adjusted for any black hole 
mass. We present all physical quantities using M = 
3.6 X IO^Mq, appropriate for Sgr A*. The data we use are 
taken from the roughly steady state portion of the simu- 
lation. We work with two different data sets. Data set 1 
consists of 32 time-slices staggered by roughly 30 hours, 
which is the orbital period at about 50 Rs- These data 
include the full 3D structure of the flow (density, pressure, 
magnetic field strength, etc.). Since much of the emission 
comes from close to the BH, however, there is significant 
evolution between these time-slices. Data set 2 thus con- 
sists of about 10,000 time-slices separated by about 30 sec, 
which corresponds to about 1/2 of the light-crossing time 
of the black hole's horizon. These data were generated 
during a run of the simulation, and only contain informa- 
tion about the synchrotron radiation not the structure of 
the flow. Because of memory constraints, we could not re- 
tain the full 3D structure of the flow at all 10,000 of these 
timcstcps. Instead, we rely on data set 1 (which has the 
full 3D structure, but sparser time sampling) to correlate 
the synchrotron emission with the structure of the flow 
and to determine optical depth effects. We use data set 2 
to probe the optically thin emission and variability on all 
relevant time-scales. 

There are several signiflcant difliculties in using the sim- 
ulations to predict synchrotron radiation. First, the den- 
sity of gas in the simulation is arbitrary because the accre- 
tion rate is not fixed in physical units; that is, although the 
simulation describes the spatial and temporal evolution 
of the density, a single normalization parameter must be 
specified to compute physical gas densities. We adjust this 
density normalization so that the total synchrotron power 
is roughly that observed from Sgr A* (this depends on 
electron temperature; see below). Given the density nor- 
malization, the total pressure and magnetic field strength 
are uniquely determined. 

The biggest uncertainty in our analysis is that the sim- 
ulations only evolve the total pressure, while synchrotron 
radiation is produced by the electrons. Since RIAFs are 
coUisionless plasmas, there is no reason to expect equal 
electron and ion temperatures or a thermal distribution of 
electrons (e.g., Quataert 2003). This is problematic be- 
cause the synchrotron emission is sensitive to the exact 
electron distribution function, about which we have no in- 
formation from MHD simulations. We account for this 
uncertainty as best as we can: we model the electrons us- 
ing a thermal distribution with a non-thermal power-law 
tail. Given a total ion -|- electron temperature (pressure) 
Ttot = Ti + Te from the simulations, we determine the 
electron temperature using T^/Ttot as a free parameter of 
our analysis. We show results for several values of Te/Ttot, 
focusing on Te/Ttot = 1 and TJTm = 1/4. This roughly 
brackets the range of electron temperatures in recent semi- 
analytic models of emission from Sgr A* (e.g., YQN03). It 
is important to stress again that there is no compelling rea- 
son to expect the electron temperature to be everywhere 
proportional to the total temperature, but this is the best 
that we can do without solving separate ion and electron 
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energy equations or carrying out kinetic simulations. 

The choice of density normahzation and electron tem- 
perature parameterization have a nontrivial impact on the 
predicted spectrum of the accretion flow. Figure 1 shows 
the optically thin synchrotron emission from a uniform 
plasma with Ttot = 2 x 10^^ K for different choices of 
Te/Ttot- For each calculation, we have adjusted the den- 
sity to fix the total synchrotron power. To understand 
these results analytically, note that the peak frequency for 
relativistic thermal synchrotron emission (where 7 oc Tg) 
scales as 



I/peafe CX oc g5 [T^/Ttotf (1) 

and the total synchrotron power scales as 

Ptotai B^j^rie (X g2 {Te/Ttotf cc q^/^iypeak (2) 

where q is the density normalization parameter (i.e., rig oc 
q) and we assume fixed /? = P/{B^/8Tr). 

Given these uncertainties, our strategy is to choose sev- 
eral parameterizations for the electron temperature (via 
Tg/Ttot) and to adjust the density normalization such that 
the total synchrotron power is reasonable. Physically, this 
corresponds to electron temperatures ~ 10^° — 10^^ K near 
the black hole and gas densities ~ 10^ — 10^ cm~^ (see also 
the one-dimensional analytic models in YQN03).^ Because 
the simulations cannot uniquely predict the total syn- 
chrotron power or peak synchrotron frequency, our analy- 
sis focuses on observables that are less sensitive to the over- 
all normalization of the emission, namely the variability 
and polarization. As we show later, the semi-quantitative 
conclusions of this paper do not depend that sensitively on 
the uncertainties in T^/Ttot and g, though detailed quanti- 
tative comparisons between simulations and observations 
are clearly premature. 

In addition to considering a thermal electron population, 
we also account for a non-thermal tail in the electron dis- 
tribution function. This power-law component is strongly 
motivated by models for the IR and X-ray emission from 
Sgr A* (e.g., Markoff et al. 2001; YQN04), and also is 
expected theoretically in models of coUisionless plasmas. 
We assume that the power-law tail has a fraction r] of the 
electron thermal energy and a power-law index p where 
n(7) CX 7~^. We fix p = 3 and choose either 77 = (no 
power-law component) or 77 = 0.05 (see YQN03). We do 
not vary p or 77 in time or space, although such variation 
is likely inevitable due to transient events such as shocks 
and reconnection (and is certainly suggested by the IR and 
X-ray flaring from Sgr A*). As a result, our calculations of 
variability at high frequencies are likely lower limits, since 
they do not include variations in the accelerated electron 
population. 

We compute the total synchrotron emission of data sets 
1 and 2 assuming ultrarclativistic electrons (following Pa- 
cholczyk 1970, henceforth P70). We have checked that 
using the transrelativistic fitting formula of Mahadevan et 
al. (1996) yields almost identical results to the ultrarcl- 
ativistic calculation for our problem. Given uncertainties 
in our viewing angle towards the accretion flow, we use a 

For the parameters chosen here, neglecting radiation in the dynamics of the accretion flow is an excellent approximation. Moreover, the 
synchrotron cooling time for the bulk of the electrons (those with 7 ~ kT^/m^c^) is usually much longer than the inflow time, so that the 
electron temperature is determined by a balance between heating and advection, rather than heating and cooling. As a result, the only way to 
model the electron energetics better than we have here is to calculate a separate electron energy equation. Were the electron cooling time very 
short (appropriate at much higher densities than we consider here), one could calculate the electron temperature everywhere in the accretion 
flow by locally balancing heating and cooling, without worrying about the Lagrangian derivative of the electron internal energy. 



simple angle averaged emissivity to compute the emission 
from data set 2 (the high time sampled data set) . We have 
assessed the validity of this approximation by computing 
the angle-dependent emission from data set 1 (which has 
the full 3D structure): the true flux is typically similar to 
that given by the angle-averaged formula to within « 50% 
(much less than the theoretical uncertainties due to the 
electron distribution function). The disagreement can be 
somewhat worse at very high frequencies {u ^ 10^'' Hz) 
when one is on the exponential tail of the synchrotron 
emissivity. 

In addition to calculating the total flux, we also calculate 

the linear polarization of synchrotron radiation in the sim- 
ulations, again assuming ultrarclativistic electrons (P70; 
eqn 3.39). We also present several calculations including 
the effects of synchrotron self-absorption. For simplicity, 
we only calculate self absorption for viewers in the equato- 
rial plane of the accretion flow. For our typical parameters 
{Tg/Ttot and q), self-absorption effects become important 
below frequencies of a few hundred GHz. We account for 
self- absorption in the thermal component by calculating 
the optical depth along a ray using the emission coeffi- 
cient and the thermal source function S'y = 2kTch''^/cP. 
The intensity along each ray is then given by adding up 
the optically thin emission from each grid point along the 
ray, weighted by exp(— r). Note that power-law electrons 
are omitted from this calculation, but each polarization 
mode is calculated separately so we can determine the po- 
larization at both optically thick and thin frequencies. 

We assume that at each time the observed emission is 
given by the steady state solution to the radiative transfer 
equation given the temperature, density, etc. of the accre- 
tion flow at that time. This is likely to be a reasonable 
approximation except for very short time-scale variability, 
and for emission from radii where the velocity of the gas 
approaches the speed of light (these limits are, of course, 
of considerable interest, but relaxing the steady state as- 
sumption is a significant effort beyond the scope of this 
work). In keeping with the Newtonian dynamics in the 
simulations used here, we also do not include gravitational 
redshift or Doppler shifts. 

It is important to point out the effects of finite spatial 
resolution on our calculations. High frequency thermal 
synchrotron emission is produced primarily at small radii 
close to the black hole, where the temperature and mag- 
netic field strength are the largest. The spatial resolution 
in this region is 0.5iis, so that at the last stable orbit 
(3i?s), there are a rather small number of grid points. We 
find that the thermal synchrotron flux at high frequencies 
(at and above the thermal peak) and the location of the 
peak are often determined by emission from a small num- 
ber of grid points. As a result, the spatial resolution is 
not quite sufficient to accurately calculate the highest fre- 
quency synchrotron emission. This is likely to be a generic 
problem in calculating emission from simulations given res- 
olution constraints. It is particularly acute above the ther- 
mal peak because the emission there is exponentially sen- 
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sitive to B and Tg. By contrast, below the thermal peak 
the emission is always well determined because it comes 
from a larger range of radii. As a check on how sensitive 
our results are to the inner few Schwarzschild radii, we 
have assessed how the emission from the flow changes if 
we only consider emission from outside (say) 4 or 6 i?^. 
We find that, although the overall amplitude of the flux at 
very high frequencies decreases, our primary conclusions 
regarding fractional variability and polarization (as enu- 
merated in §4) are essentially unchanged. Thus we believe 
that these results accurately reflect the underlying dynam- 
ics of the accretion flow. 

3. RESULTS 

To begin we present results for unpolarized synchrotron 
emission, witliout taking into account the effects of syn- 
chrotron self-absorption. We then discuss the polariza- 
tion of the emission and the effects of synchrotron self- 
absorption. For some of the frequencies considered below 
it is not self-consistent to ignore self-absorption, but we do 
so initially to illustrate general points about synchrotron 
emission from the flow. 

3.1. Optically Thin Emission 

To become familiar with the emission it is useful to ex- 
amine a time-slice and 'see' what the flow would look like 
if it could be resolved with an array in the sub-mm. This 
is shown in Figure 2 for a typical time-slice, as viewed 
from an arbitrary angle off of the plane at « 450 GHz, 
which is near the peak in the thermal emission (see Fig. 
7). We only show the central 16 which is by far the 
dominant emission region (the intensity plot in Fig. 2 is 
logarithmic). Most of the emission originates close to the 
central hole, in the equatorial plane, rather than, for in- 
stance, in polar outflows or jets, which do not appear in 
this simulation. This structure changes somewhat from 
time-slice to time-slice but the prominence of an equato- 
rial region near the hole is common to all. At lower fre- 
quencies, where v ^ Vpeak, large radii are more important 
and so the emission is much less centrally concentrated. 
Note that the caleiilations shown in Figure 2 do not in- 
clude any General Relativistic photon transport (or even 
obscuration by the BH), which may give rise to unique 
signatures in the observed image (e.g. Falcke, Melia, & 
Agol 2000). Nonetheless, these calculations are encourag- 
ing because the emission arises so close to the hole, where 
GR effects would indeed be important. 

Figure 3 shows the variability of the emission with time 
at three representative frequencies, from the radio to the 
IR. The calculations are for Te/Ttot = 1/4 and a power- 
law tail with -q = 0.05. The = 10 GHz and = 80 GHz 
emission are produced by the thermal component, while 
the = 81.9 THz emission (K-band IR) is produced by 
the non-thermal tail. The most striking feature of these 
results is the strong variability of the emission with time. 
The flux varies by up to an order of magnitude, often on 
time-scales as short as an hour. Figure 3 also shows that 
the amplitude of the variability increases, and the time- 
scale for variability decreases, at higher photon frequen- 
cies. This is because higher photon frequency emission 
originates closer to the hole where the dynamical time- 
scales are shorter. In addition, we flnd that variability of 



the flow parameters in the inner region is larger than the 
variability at large radii, accounting for the fact that the 
higher frequency emission that originates close to the hole 
is more variable. 

To demonstrate the origin of this variability. Figure 4 
shows how the flux at 100 GHz varies on long time-scales, 
as compared with the magnetic fleld strength {\B\), tem- 
perature (Te), and density (rie), averaged from 2 to 6 iJg. 
It is clear that the variability is primarily due to changes in 
the magnetic field strength, with fluctuations in temper- 
ature and density being less important (the correlation 
coefiicients between vLi, and B, n, and T, are 0.93, 0.66, 
and 0.57, respectively). Note that because the flow has 
/3 ~ 10 — 100, even large changes in the magnetic field 
strength do not necessarily lead to appreciable pressure or 
density changes, though they do cause significant variabil- 
ity in the synchrotron emission. 

In the previous section we emphasized that the predicted 
emission depends on how we parameterize the electron 
temperature and distribution function. Figure 5 contrasts 
these different parameterizations in terms of the variability 
of the flow on different time-scales and at different photon 
frequencies. Figure 5 displays the power in Fourier com- 
ponents corresponding to one day and one hour periods, 
normalized to the mean flux (see caption for details). On 
hour time-scales the two models with purely thermal elec- 
trons look very similar, except that the variability at a 
given frequency is smaller for the model with hottc;r elec- 
trons. The reason for this is that higher electron tempera- 
tures imply that at a given photon frequency the emission 
arises from larger radii (so that B is smaller), where the 
dynamical timescales are longer and thus the short time- 
scale (~ hour) variability is weaker. 

Figure 5 also shows that higher photon frequency emis- 
sion is significantly more variable than lower photon fre- 
quency emission on hour time-scales. This is because the 
lower photon frequency emission includes significant con- 
tributions from large radii where the orbital time-scales 
that govern the variability of the flow are » 1 hour. The 
trend of increasing variability with increasing photon fre- 
quency is present, but noticeably less dramatic, on day 
time-scales. This is simply because most of the radii that 
contribute to the observed emission have orbital periods 
^ 1 day, and thus are significantly variable on ~ day time- 
scales. Finally, Figure 5 shows that non-thermal emission 
begins to dominate thermal emission at v « 10^^ — lO^'"* 
Hz. In addition, the models with pure thermal emission 
are more variable than those with non-thermal emission 
at these frequencies: this is because ^ lO'^^ — 10"'^^ Hz is 
near or above the peak in the thermal synchrotron emis- 
sion (see Fig. 7), and so the flux is exponentially sensitive 
to variations in the flow parameters. Of course, the flux 
(and the variable flux) in the non-thermal tail at high fre- 
quencies is much larger than that produced by the thermal 
electrons, but the normalized variability is smaller. 

3.2. Linear Polarization 

We have also computed the net linear polarization of 
synchrotron emission in the simulation (neglecting Fara- 
day rotation and optical depth effects; the latter are con- 
sidered in the next section). The results are shown as a 
function of frequency in Figure 6 for all of the time-slices in 
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data set 1 (each of which is separated by about 30 hours) ; 
for these calculations we chose to view the flow 30 degrees 
off of the equatorial plane, set T^jTtot = 1/4, and included 
a non-thermal tail with t] = 0.05. As Figure 6 shows, the 
magnitude of the linear polarization is substantial, typi- 
cally 30 — 40%. There is some variation with frequency, 
though it is not particularly dramatic. The polarization 
vector lies perpendicular to the equatorial plane. It arises 
due to the coherent toroidal magnetic field in the flow, 
analogous to the polarization of synchrotron radiation in 
the Galaxy (e.g.. Beck 2001). The polarization would be 
larger {^i 60%; see Fig. 8) if we viewed the flow edge- 
on, while it would vanish for face-on viewing angles. Note 
also that there is appreciable variability in the magnitude 
of the polarization, roughly order unity changes on day to 
week time-scales. The variability is somewhat larger at 
high frequencies, again because this emission arises closer 
to the black hole. 

3.3. Synchrotron Self-Absorption 

In the previous sections we have assumed that the syn- 
chrotron emission is optically thin at all frequencies of in- 
terest. This assumption is incorrect at low frequencies 
where synchrotron self- absorption becomes important and 
the emission becomes optically thick. In this section we 
present initial results on the emission and variability at 
optically thick frequencies, leaving a more comprehensive 
discussion to future work. We have computed the total flux 
as a function of frequency including self-absorption for the 
case when the flow is viewed in the equatorial plane, per- 
pendicular to the rotation axis (sec §2 for the method); we 
considered thermal electrons only and took Tg/Ttot — 1/4. 
The top panel in Figure 7 shows the average and RMS 
emission both with (solid) and without (dashed) optical 
depth effects, while the bottom panel shows the normal- 
ized RMS variability for these two cases. ^ The averages 
and RMS are taken over all 32 time-slices (data set 1) 
for which we have the full 3D structure necessary to com- 
pute optically thick emission. Since these time-slices are 
separated by ~ 1 day, we do not have good constraints 
on short time-scale variability for the optically thick emis- 
sion. Note also that the RMS variability plotted in Figure 
7 yields a quantitative measure of variability that is sys- 
tematically smaller (by a factor of few) than the Fourier 
analysis shown in Figure 5 (compare with the f = 1 / day 
plot), despite the fact that the data are essentially the 
same. The Fourier amplitude is a better measure of the 
'by-eye' variability, but with only 32 (somewhat unevenly 
distributed) time-slices a Fourier-amplitude version of Fig- 
ure 7 was not feasible. 

Figure 7 shows that self-absorption becomes important 
at ^ 10^^ — 10^^ Hz and that the emission is, of course, 
substantially suppressed below this frequency. Interest- 
ingly, Figure 7 also shows that the fractional variability is 
also suppressed below the self-absorption frequency. The 
reason is simple: optically thin calculations erroneously in- 
clude emission from small radii, where the flow is strongly 
variable. By contrast, when optical depth effects are in- 
cluded the photosphere moves out to ~ 10 — 50i?s (depend- 

^ To check whether the optically thick results at low frequency are sensitive to spatial resolution (e.g., resolving the r ~ 1 surface), we 
interpolated the flow structure to a finer grid and recalculated the emission. The results were essentially identical to those in Figure 7 for 
V ^ 10^^ — lO'^^ Hz. The higher frequency (optically thin) emission did change somewhat, however, because it is dominated by very small radii 
(see the discussion at the end of §2). 



ing on frequency). The net effect is that the fractional 
variability is smaller when synchrotron self-absorption is 
accounted for because the emission from small radii is ex- 
cluded. 

We have also calculated the polarization of the syn- 
chrotron emission including optical depth effects (Fig. 8). 
We neglect the eff'ects of Faraday rotation in this cal- 
culation and again consider only thermal electrons with 
Te/Ttot = 1/4. These calculations, in contrast to the op- 
tically thin polarization calculations shown in Figure 6, 
place the observer in the equatorial plane of the disk, thus 
maximizing the observed polarization. As Figure 8 shows, 
the dominant effect of self- absorption is to substantially 
decrease the polarization fraction at self-absorbed frequen- 
cies ^ 10^^ Hz. This is because emission from the central 
optically thick impact parameters dominates the emission 
from peripheral optically thin impact parameters, and op- 
tically thick thermal emission is unpolarized. Figure 8 also 
shows that the variability of the polarization fraction re- 
mains largely unchanged when optical depth effects are 
taken into account. The angle of the polarization is also 
unchanged, since it is still determined by the predomi- 
nantly toroidal field. 



4. DISCUSSION 

Our results on synchrotron radiation from MHD simu- 
lations of radiatively inefficient accretion fiows can be suc- 
cinctly summarized as follows: (1) the emission is highly 
variable: the fiux can change by up to an order of magni- 
tude on time-scales as short as the orbital period near the 
last stable orbit (~ 1 hour for Sgr A*). (2) the variabil- 
ity is stronger and more rapid at high photon frequencies 
because the high frequency emission arises from closer to 
the black hole; (3) The emission at self-absorbed frequen- 
cies is less variable than the emission at optically thin 
frequencies. This is in part a consequence of point (2); 
optically thick emission arises further from the hole where 
the variability is weaker and less rapid. (4) At optically 
thin frequencies, the synchrotron emission is linearly po- 
larized at the ~ 20 — 50% level (unless the fiow is viewed 
face-on or there is significant Faraday depolarization); the 
polarization vector is perpendicular to the equatorial plane 
of the accretion flow and is due to the coherent toroidal 
magnetic field. The magnitude of the polarization is it- 
self significantly variable (by factors of few). (5) For a 
roughly thermal electron distribution function, the polar- 
ization fraction decreases significantly at optically thick 
frequencies. 

We believe that these results are reasonably robust in 
spite of several significant uncertainties in our analysis. In 
particular, how we treat the electron temperature in the 
fiow does not significantly change these conclusions (see 
Fig. 5). We have also checked that our conclusions arc 
unchanged if we only consider emission from gas outside 
of (say) 4i?s (neglecting emission from between 2 — 4 iig). 
This is encouraging because our pseudo-Newtonian treat- 
ment of the dynamics is particularly inaccurate at very 
small radii. We suspect that the results presented here 



6 



are characteristic of models in which the emission is domi- 
nated by gas in the accretion flow itself. If, however, there 
is preferential electron acceleration/heating in a low den- 
sity corona or in an outflow, this could lead to variabil- 
ity and polarization signatures significantly different from 
those presented here. A much more careful treatment of 
the electron thermodynamics (beyond MHD) is required 
to assess this uncertainty. In future work, it would also 
be interesting to carry out the radiative transfer including 
Faraday rotation and General Relativistic photon trans- 
port. 

4.1. Application to Sgr A* 

Wc suggest that the conclusions highlighted above re- 
garding synchrotron emission in the MHD simulations are 
reasonably consistent with observations of Sgr A* from the 
radio to the X-ray (summarized in §1). For example, our 
calculations produce significant linear polarization in the 
mm-IR emission, in accord with observations (Bower et al. 
2003; Genzel et al. 2003). We also find a strong decrease in 
the linear polarization fraction with frequency below ~ 100 
GHz due to optical depth effects. It is probably not sur- 
prising that the precise values of the observed polarization 
fractions arc not recovered in our calculations; we do not as 
yet include Faraday depolarization, which should further 
decrease the polarization fraction (Quataert & Gruzinov 
2000), nor can we easily constrain the overall density of 
the flow, which determines the exact frequency at which 
optical depth effects set in. 

The calculations presented here predict that the polar- 
ization angle is perpendicular to the equatorial plane of the 
accretion flow, and should not change signiflcantly with 
frequency. The former prediction is somewhat tricky to 
test, however, because it is not obvious what the orien- 
tation of the accretion flow is on the sky. The Bardeen- 
Petterson effect may not be efficient for thick accretion 
disks (e.g., Natarayan & Pringle 1998), so it is unclear 
whether the angular momentum of the flow at small radii 
is tied to that of the hole. If not, the orientation of the flow 
may be set by the angular momentum of the stars whose 
winds feed the hole (see Levin & Beloborodov 2003 and 
Genzel et al. 2003a for discussions of the stellar angular 
momentum). If, however, the disk's angular momentum is 
tied to the hole's, then the orientation of the flow at small 
radii will depend on the angular momentum accreted by 
the hole over its lifetime, which is not known. 

Our calculations naturally produce factors of ~ 10 vari- 
ability in the IR emission on ~ hour time-scales (Fig. 3), 
and thus may account for some of the observed IR vari- 
ability from Sgr A*. They do not, however, quite produce 
siifficient variability on time-scales as short as is observed, 
^ 10s of minutes (less than the orbital period at the last 
stable orbit for a non- rotating black hole). This could be 
because our calculations have not been carried out for ro- 
tating black holes or because the most rapid IR variability 
traces particle acceleration, not the accretion flow dynam- 
ics.^ It is worth noting that we do not see any evidence for 
quasi-periodic oscillations in our calculations, though we 
have not included Doppler boosting that could modulate 
the emission on the orbital time-scale (e.g., Melia et al. 

^ Note also that our spatial resolution is the poorest at the small radii 
(e.g., turbulent eddies) that vary on ~ 10 minute time-scales. 



2001). 

More generally, our calculations produce variability that 
increases in amplitude and decreases in time-scale with in- 
creasing frequency (Figs. 3 & 5), as is also observed from 
Sgr A*. In particular, the suppression of variability at 
optically thick frequencies (Fig. 7) may contribute signifl- 
cantly to the differences in the observed variability in the 
radio, mm, and IR (it is also possible that this is due to 
a change in the dynamical component responsible for the 
emission; e.g., a jet becoming important at lower frequen- 
cies as in Yuan et al. 2002). In addition, although we have 
not calculated X-ray emission, the large amplitude, hour 
time-scale variability we see from gas close to the black 
hole is somewhat reminiscent of the X-ray flares observed 
by Chandra and XMM. Moreover, since some of the X-ray 
flares could be produced by synchrotron emission (e.g., 
YQN04), it is plausible to suppose that the synchrotron 
variability we calculate here could extend to higher fre- 
quencies as well. Finally, it is worth reiterating that that 
the calculations presented here are also likely lower lim- 
its to the variability ai v ^ 10^^ Hz because we do not 
include variations in the non-thermal electron population 
or effects from synchrotron sclf-Compton emission. Such 
transient particle acceleration appears to be required to 
explain the very large amplitude X-ray flares (Markoff et 
al. 2001; YQN04). 

Our calculations predict that variability at different fre- 
quencies should be strongly correlated so long as both fre- 
quencies are optically thin (see Fig. 3). We also pre- 
dict that the time delay between the emission at differ- 
ent optically thin frequencies should be quite small, ^ 
hours. These predictions may be testable by correlating 
sub-mm and IR variability from Sgr A*. Note, however, 
that strong temporal or spatial variations in electron ac- 
celeration could modify these predictions. Our results also 
suggest that emission at optically thick frequencies should 
not be as well correlated (though we have not shown this 
explicitly). The reason is that different frequencies then 
probe different radii; since the turbulence at one radius 
need not be well correlated with the turbulence at another 
radius, the same follows for variations in the synchrotron 
emission. 

In spite of a few shortcomings, the general agreement 
between the properties of synchrotron emission in our cal- 
culations and observations of the Galactic Center is en- 
couraging. It supports a model in which much of the high 
frequency emission from Sgr A* is generated by a turbu- 
lent magnetized accretion ffow close to the black hole, and 
encourages more reflncd calculations of emission from nu- 
merical simulations of RIAFs. 
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re the IR emission is produced, so we may not fully resolve structures 
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Fig. 1. — uLi, from an arbitrary volume of gas with rie = g X 10® cm^'', B = y/q X 70 G. and Te = 2 X 10^^ {Tc/Ttot) K, for various values 
of Tfi/Ttot and g; note that B oc y'g is equivalent to /3 = const. In each case the density parameter q has been adjusted so that the total 
synchrotron power is roughly the same: q Ri 2.49, 6.23, 7.56, 10, & 12.5 for Te/Ttot = 1, 2/5, 1/3, 1/4, & 1/5. 




Fig. 2. — An image of the central 16 Rb at a particular time-slice for u 450 GHz, near the peak of the thermal emission. The intensity 
scale is logarithmic and the brightest area is in the equatorial plane near the last stable orbit, where the density, temperature, and magnetic 
field strength are the largest. For u <^ 450 GHz (below the thermal peak), the emission is much less centrally concentrated. 
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Fig. 3. — Optically thin light-curves at three representative frequencies, for T^/Ttot = 1/4 and ?? = 0.05, where r] is the fraction of the 
electron thermal energy in a power-law tail with n(7) oc 7~^. The emission at v = W and 80 GHz is produced by the thermal electrons, while 
the V = 81.9 THz emission (K-band IR) is produced by the power-law electrons. 
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Fig. 5. — EELch panel displays variability for 3 different models of the electron thermodynamics as a function of photon frequency {rj is the 
fraction of the electron energy in a power-law tail) . The top panel shows hour time-scale variability while the bottom is for day time-scale 
variability. The Y-axis is a measure of the normalized variability on a particular time-scale, where (P) is the time average of fL^, and 
P (/) = \ {uLv{t)) I, where is the Fourier transform. These calculations assume optically thin synchrotron emission; optical depth effects 
are considered in Fig. 7. 
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Fig. 6. — Optically thin lineaj: polarization fraction as a function of frequency for each of the 32 time-slices for which we have the full 3D 
structure (see Fig. 8 for optical depth effects). Each time-slice is separated by « 30 hours. All curves assume Te/Ttot = 1/4 and rj = 0.05. 
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Fig. 7. — The top panel displays the mean {vLv) and the RMS variability of the emission including (solid line) and excluding (dashed line) 
optical depth effects. The bottom panel shows the variability as a fraction of the mean. Note that the RMS shown here yields a quantitative 
measure of variability that is systematically smaller (by a fax;tor of few) than the Fourier analysis shown in Figure 5 (see text for more 
discussion) . 
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Fig. 8. — The top panel shows the average polarization fraction (11) of the emission including (solid line) and excluding (dashed line) 
optical depth effects (averaged over all 32 time-slices for which we have spatial data). The bottom panel shows the RMS variability of the 
polarization fraction. These calculations assume purely thermal electrons with T^ITtot = 1/4 and that the disk is observed edge on. 



